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ABSTRACT 



We present evidence to support an earlier indication that the Galaxy is embedded in an 
extended, highly inclined, triaxial halo outlined by the spatial distribution of companion galaxies 
to the Milky Way. Signatures of this spatial distribution are seen in 1) the angular variation of 
£^ | the radial-velocity dispersion of the companion galaxies, 2) the spatial distribution of the M 31 

■ sub-group of galaxies, 3) the spatial distribution of the isolated, mainly dwarf irregular, galaxies 

of the Local Group, 4) the velocity anisotropy quadrupole of a sub-group of high- velocity clouds, 
and 5) the spatial distribution of galaxies in the Coma-Sculptor cloud. Tidal effects of M 31 
and surrounding galaxies on the Galaxy are not strong enough to have affected the observed 
structure. We conclude that this distribution is a reflection of initial conditions. A simple 
galaxy formation scenario is proposed which ties together the results found here with those of 
Holmberg (1969) and Zaritsky et al. (1997) on the peculiar distribution of satellites around a 
| large sample of spiral galaxies. 
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52 ' 1< Introduction 

C3 ■ 

Our present view of the Milky Way had its origin over 80 years ago with the classic work of Shapley 
(1918). The 'Shapley Galaxy' had the Sun offset from the center in a disk embedded in a roughly spherical 
halo of globular clusters. At the time of Shapley 's work the only available probes of the outer halo were 
globular clusters. It was not until 20 years later that the first dwarf spheroidal companion to the Galaxy 
was discovered, also by Shapley. A large increase in the discovery of distant low surface brightness satellites 
came from plates taken with the 48-inch Palomar Schmidt telescope in the early 1950's and has continued 
up to the present time. It was then not until the mid 1970's that almost simultaneously Lynden-Bell (1976) 
and Kunkcl & Demers (1976) noted that many of the outlying satellites lay in a plane nearly coincident with 
the Magellanic Stream of neutral hydrogen. Since then other streams have been discussed (c.f. Lynden-Bell 
1982, Majewski 1994, Lynden-Bell & Lynden-Bell 1995). In what follows we consider this problem from 
a different perspective, clarify and extend our previous work, (Hartwick 1996a, b, hereafter H96a,b), and 
reach similar conclusions using an updated sample and an improved method of data analysis. There it was 
argued that our 'Shapley Galaxy' is embedded in an elongated, highly inclined outer halo defined by the 
spatial distribution of the associated companion galaxies to the Milky Way. We strengthen this picture by 
showing that a very similar spatial distribution applies to the M 31 sub-group of galaxies and to the Local 
Group as defined by the outer isolated, predominately dwarf irregular, galaxies. In addition, we conclude, 
on the basis of an apparent similarity in these spatial distributions with nearby larger-scale structure, that 
the results may be understood as a reflection of initial conditions. 
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2. Defining the Outer Halo of the Galaxy 



2.1. Sample Selection 



Any method of defining the limits of the outer halo will be somewhat arbitrary and for this reason we 
look for common properties within different subsamples of both nearest companion galaxies and globular 
clusters with Galactocentric distances between 25 kpc and 450 kpc. Within these limits are 10 companion 
galaxies (treating the LMC and SMC as one system) and 15 globular clusters. The data on the following 
satellite galaxies were taken from the recent review by Mateo (1998): LMC, Fornax, Leo I, Sculptor, 
Sextans, Leo II, Ursa Minor, Carina, Draco and Phoenix. The globular cluster data were taken from the 
most recent web-based update of the compilation of Harris (1996) and included the following clusters: 
AM-1, Eridanus, Pal 2, NGC2419, Pyxis, Pal 3, Pal 4, AM-4, NGC5694, NGC5824, Pal 14, NGC6229, Pal 
15, NGC7006, and Pal 13. No radial velocities were available for Pyxis or AM-4 and an 'average' value of 
16 km sec -1 was assumed for Phoenix. 



Whereas in the previous work (H96a,b) the spatial distribution was determined by fitting an ellipsoid 
to the data by least squares, in this work we have characterized the spatial structure by a symmetric tensor 
with average values XiXi/di,yiyi/di,ZiZi/di,Xiyi/di,XiZi/di,yiZi/di as components, where di is the distance 
of the i th object from the origin (in this case the Galactic center with Ro = 8 kpc assumed) and Xi,yi, 
and Zi are the projections of di along the principal axes of the Galactic coordinate system. Eigenvalues of 
the matrix were obtained in the standard way and are given below along with their directions for various 
subsamples of the data. Throughout this paper the quoted one sigma errors were computed by the bootstrap 
method. The errors on quantities involving a direction (i.e. the spatial and radial velocity eigenvalues), 
were computed in two ways. The first allowed for errors due to differences in direction (as below) while the 
second assumed a fixed direction and projected the resulting dispersions onto the principal axes. The errors 
computed by the second method were found to be so similar to those from the first that we have not given 
them here. The results of these procedures are given below where the units of the eigenvalues (e a ,6, c ) are 
kpc, and I and b are Galactic coordinates expressed in degrees. 



2.2. The Spatial Distribution 
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Galaxies — Phoenix & Leo I (n=8) 
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The above solutions indicate that the outer satellites define a triaxial spatial distribution whose highly 
inclined long axis is nearly independent of the outer most object in the sample. These directions differ 
by only 18° ± 14 [(l)-(2)] and 8.2° ± 12 [(l)-(3)]. The directions of the other principal axes are also 
quite similar. In spite of the fact that the number of objects outlining this structure is small, our working 
hypothesis is that the above result defines a real structure and in what follows we look for supporting 
evidence. We start by considering the spatial distribution of the 15 outer-most globular clusters and find: 

Globular Clusters (n=15) 
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The above result suggests that there may be a real difference between the spatial distributions of the 
clusters and the galaxies. The cluster distribution appears to be more nearly oblate and not as flattened 
(e c /e a = 0.46^q'J|) as the galaxy distribution. The inner cutoff of ~ 25 kpc in cluster selection is not a 
significant factor as very similar results are obtained if we omit the radial cutoff and instead choose clusters 
by their [Fe/H] values. For 99 clusters with [Fe/H] < — 1.0 the shortest axis of the still nearly oblate spatial 
distribution is directed towards I — 307^111, b = —3.31^15.1 with axis ratio e c /e a = 0.55io'n- 

While both of the above samples are discouraging small, there is the strong suggestion that both 
luminous components of the outer halo of the Galaxy form structures whose minor axes are highly inclined 
to the Galactic rotation axis. This result is similar to that found in the seminal work of the mid 1970's. For 
comparison with this earlier work, the pole of the Magellanic Stream lies at I = 11°, b = —10°, (Lyndcn-Bcll 
1976) and the pole of the 12 members of the 'Magellanic Plane Group' of Kunkel & Demers (1976) is at 
I = 347°, b= 22°. 

As a further check on our procedure we have examined the spatial distribution of the 34 metal-rich 
globular clusters with [Fc/H] > —0.7. The somewhat surprising results are: e a = 2.3llg'53g towards 
I = 347tl?;jj, b = — 11.0±|"J§; e b = 1.13±g;||| towards I = 77.7111;?, b = -1-891?;^ and e c = 0.343l°;J^ 4 
towards I = 357lgf ;f , b = 78.8l7;;j!g. A full discussion of these results, including observational selection 
effects, is beyond the scope of this paper, but we note that the short axis of the triaxial distribution now 
points toward high Galactic latitude while the long axis is only ~ 17° off the Sun-center line which is 
very reminiscent of the bar-like feature believed to lie towards the Galactic center (c.f. §10.2 of Binney & 
Merrifield 1998). 

Before proceeding it is reasonable to consider how configurations such as we have seen in solutions 
(1),(2), and (3) could arise in practice. The possibility that the peculiar configuration observed has arisen 
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by chance from what is actually a spherical distribution in the presence of the zone of avoidance was shown 
by a simple Monte Carlo simulation to have low probability (H96a). Two other possibilities are that we 
are seeing either the effects of dynamical evolution or a reflection of initial conditions. It is well known 
that dynamical friction will tend to bring satellites which are initially inclined to the plane of a galaxy into 
the plane working most efficiently on those objects closest to the plane initially (c.f. Quinn & Goodman 
1986). While this mechanism undoubtedly does operate in the inner part of the Galaxy where orbital times 
are relatively short, it is more difficult to understand how it could significantly alter the distribution of 
the outer-most objects (Quinn & Goodman 1986, Zaritsky et al. 1997). For now, we will consider the 
possibility that we may be seeing the effects of initial conditions, and after the following analysis of the 
radial velocities of the above samples will look for similar configurations in other nearby systems. 



2.3. The Angular Variation of the Radial- Velocity Dispersion 

In order to look for possible kinematic signatures of the peculiar spatial distribution, radial-velocity 
dispersion tensors were computed in H96b for different samples of halo objects by the method of maximum 
likelihood. Here we confine our attention to the distant samples where the observed radial velocities are 
primarily a measure of the radial component of the velocity dispersion, and we use a mathematically 
simpler and equivalent procedure analogous to that used in determining the spatial distributions. First, 
the individual heliocentric radial velocities were corrected for the solar motion with respect to the local 
standard of rest, LSR, (using the value from Binney & Merrifield 1998) and the rotation of the LSR with 
respect to the center of the Galaxy (assuming V rot = 220 km sec -1 ). These radial velocities were then 
referred to the Galactic center by multiplying each by the secant of the angle at the object between the 
Galactic center and the Sun. The motion of the Galaxy with respect to the halo objects, Vq, was then 
determined by least squares. We define a radial-velocity dispersion tensor by computing the coefficients 
V X V X , VyVy, V Z V Z , V x V y , V X V Z , V y V z where V x = V r t c x x etc and V r ^G is the observed radial velocity corrected 
for the two components of solar motion, referred to the Galaxy center and corrected for the above values 
of Vq. The direction cosines x, y, and z are once again calculated with respect to the Galactic center. As 
the distance to the Galactic center is small compared with the Galactocentric distances of our objects, we 
do not expect much sensitivity to any tangential component of the velocity dispersion. We confirmed this 
by performing heliocentric solutions which give a different weighting to any tangential component and as 
expected found results almost identical to those below. Both methods yield large but similar values for Vq- 
As the direction of Vq is usually close to being 90° away from the major axis we are probably seeing the 
effects of rotation but given the extremely poor sensitivity to tangential motions, and the very small sample 
numbers such a motion will not be well determined. The radial-velocity dispersion tensor calculated as 
above was then diagonalized and the results are given below along with the corresponding value of Vq in 
km sec -1 . The units of the eigenvalues denoted by a are km sec -1 . 

All Galaxies (n=10) 
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Galaxies — Phoenix (n=9) 
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Globular Clusters (n=15) 
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As can be seen from above, all three galaxy solutions are nearly perfectly aligned with the corresponding 
spatial solutions of the previous section and in the sense that the largest dispersion axis coincides with 
the longest spatial axis and vice versa. We cautiously interpret this result as support for our hypothesis. 
That this is not the case for the outer clusters, however, is another indication that the cluster component 
is different from the galaxy component. It should emphasized that what is being measured above is the 
radial- velocity dispersion and that without knowledge of the tangential component of the velocity dispersion 
the kinematical state of the outer halo remains incomplete. 



3. The Spatial Distribution of Satellites in the M 31 Sub-Group 

If our hypothesis is correct then we might expect to see similar structures around other galaxies. M 31 
would appear to be an ideal test case as it has at least 11 satellite galaxies with reasonably well determined 
distances according to the recent review of Mateo (1998). The galaxy EGB0427+63 was not included in our 
sample due to the large quoted uncertainty in its distance and an updated distance of 660 kpc was assumed 
for IC10 (Sakai et al. 1999). The spatial distribution of the 11 members of this grouping of galaxies as 
denoted by Mateo(1998) (also see Karachcntsev, 1996) was determined below using identical procedures as 
in §2.2 except that d{ is now taken to be the true distance from M 31. As before, the units of e a j,,c are kpc. 



M 31 Sub-group (n=ll) 
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Noting that the distance errors for the M 31 companions are generally larger than those for the Galaxy 
and that the low Galactic latitude of M 31 almost certainly contributes to the incompleteness of the sample, 
the results in (9) are quite remarkable not only because they suggest a similar triaxial structure around 
M 31 but that the long axis of this structure is within 16° ± 17 of that found for the Galactic satellites in 
solution (1). A solution combining both the 10 Galaxy companions and the M 31 satellites translated to a 
common origin is given below. 



Galaxy & M 31 Satellites Combined (n=21) 
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Figure 1 shows a projection of the 3-dimcnsional representation of this combined sample. As an aside 
one could also consider Figure 1 to be a visual illustration of the fact that the ratio of the density of 
luminous matter to the total matter density i.e. fl^/flo ~ 0.01 when one realizes that the pile of lighted 
baryons which we call the Galaxy is contained well within the elliptical contour at the origin. 



Fig. 1 here 



4. A Relationship with Larger Scale Structure 

Given the results of the previous sections we are encouraged to trace the signature of solutions (1), 
(2), and (3) to larger scales. We start by considering the spatial distribution of the isolated members of 
the Local Group, then consider the kinematics of the high-velocity clouds which are potential Local Group 
interlopers and finally look at the distribution of galaxies in the Coma-Sculptor cloud, a nearby large-scale 
structure in which the Local Group is embedded. 



4.1. The Spatial distribution of Local Group Galaxies 

Excellent reviews of the properties of Local Group galaxies have recently been given by Mateo (1998) 
and van den Bergh (2000). For our purposes we focus attention on the spatial distribution of a subset of 
these galaxies and consider only those galaxies which are relatively isolated (principally the dwarf irregulars, 
i.e. those galaxies which are not obvious satellites of/nor either the Galaxy or M 31). Given that this 
gas-rich subset likely contains potential but as yet 'unincorporated' galaxy building blocks, we argue that 
there is the reasonable possibility that the spatial distribution of these galaxies could provide us with an 
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outline of the original shape of the Local Group. Mateo (1998) argues the case for a more extended Local 
Group membership than van den Bergh, and based on the discussion in Mateo (1998) we have chosen 
our sample from his compilation. It includes: WLM, NGC55, IC1613, Leo A, NGC3109, GR8, SagDIG, 
NGC6822, DDO210, IC5152, Tucana, UKS2323-326, and Pegasus. Mateo's (1998) distances were used for 
each of the objects. The center of mass was assumed to be at 460 kpc from the Galaxy in the direction of 
M 31. The procedure followed is identical to that used to determine the spatial distribution of Galactic 
satellites, and the results are given below. The units of e a .b, c are kpc. 

Local Group (Isolated Dwarfs) (n=13) 
e a = 654ii« i = 349±?§;| 6 = -54.7™ 

e b = 335±^ 5 1 = 307±U:f 6 = 28.0™ (H) 
e c = 128™ ^48.2^o.i 6=19.7™ 

In spite of the larger uncertainties in the distances of many Local Group galaxies, we note the similarity of 
the above distribution to solutions (1), (2), and (3) for the Galaxy (except for a reversal of the intermediate 
and smallest axes). Had we chosen our subset from van den Bergh's (2000) compilation, NGC55, NGC3109, 
GR8, IC5152, and UKS2323-326 would be deleted from the list above. A solution with the remaining 
8 galaxies also shows a triaxial distribution with the long axis at I — 358^1^, b = — 34.ll20.9- Due 
to uncertainties both in defining the zero-velocity surface, which leads to the definition of Local Group 
membership, and in defining the 'isolated, gas-rich' subset which by necessity involves small numbers, we 
are prompted to look beyond the Local Group and consider the spatial distribution of the galaxies making 
up the nearest large-scale structure. Before doing so, we examine the kinematics of the high-velocity clouds. 

4.2. Kinematics of High- Velocity Clouds 

It is clear that additional probes of the region occupied by the satellite galaxies are desirable in order 
to strengthen the hypothesis. Recent work on high-velocity clouds (HVC's) attempts to make the case 
for them being Local Group interlopers (c.f. Blitz et al. 1999, Braun & Burton 1999). Solar motion 
solutions by both groups of authors reveal a residual which is interpreted as a systemic infall towards the 
Sun (c.f. Figure 5 of Braun & Burton (1999)). Below we examine the angular variation in the velocity 
field by evaluating the quadrupolar term using the isolated, compact (but incomplete) sample of Braun & 
Burton. The starting point is the following expression for the radial-velocity component of the velocity field 
expanded about the origin (the Sun in this case) for an individual HVC 

v r = n • vo + n • D • r (12) 

where n = r/\r\ is the unit vector containing the object's direction cosines, — vq is the solar motion, D is 
the displacement tensor (c.f. eqn 2.72 of Ogorodnikov, 1965) and r is the position vector. Since only the 
directions and not the distances are known for the clouds, we replace r with n in the above equation and 
solve for the components of v and the symmetrical anisotropy tensor by least squares. Eigenvalues are in 
km sec -1 . 
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Velocity Anisotropy Tensor for HVC's (n=66) 
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The monopole term is ^„ ono = ~105lg' 3 Q km sec -1 and this value (the average of the original 3 
eigenvalues) has been subtracted from the original 3 eigenvalues to give the values tabulated above. Once 
again one can identify the directions in the above solution with those of (1), (2), and (3). This identification 
should be considered tentative until the sample is complete but it should also provide even more incentive 
to determine the true nature (at least the distances) of the HVC's. 



4.3. Nearby Large-Scale Structure 



We note that the shortest axis of the Local Group spatial distribution lies within 13° ± 10 of the 
direction of the supergalactic pole at / = 47.4° and b = +6.32°, suggesting a possible connection with 
nearby large-scale structure. Examination of the Nearby Galaxies Atlas (Tully & Fisher 1987) shows that 
the Local Group is embedded in a filament which the above authors refer to as the Coma-Sculptor Cloud. 
The spatial distribution of the galaxies in the Coma-Sculptor Cloud was quantified in the same manner 
described above. A solution was carried out for those galaxies with a group prefix designation of 14 from 
Tully's (1988) catalog and is given below. The center of gravity of the 214 galaxies was determined by 
giving each galaxy unit weight and was found to lie at 2.8 h^Mpc towards 1 = 155°, b = 83.5°. The units 
of e a ^, c below are h _1 Mpc with h = H /100. 



Coma-Sculptor Cloud (n=214) 
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Once again we note that the long axis of the distribution is pointing at high Galactic latitude while the 
short axis differs from the position of the supergalactic pole by only 11°. 



4.4. Tidal Effects 



Since gravity is responsible for the growth of structure in the universe it is relevant to consider the 
effects of the quadrupolar (tidal) component of the gravitational field on the structures we have discussed 
above. Following Raychaudhury & Lynden-Bell (1989) (RLB89), we define the tidal tensor as 
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where the tensor components £1,1 = 3 x Ji Xa — l,ii,2 = 3 x ^ Va , f 1)3 = 3 ^g z ° etc and as before d a is the 
distance of the a th object from the origin under consideration and x a ,y a and z a are the components 
of d a projected onto the principal axes of the Galactic coordinate system. M a is the mass of the a th 
galaxy. In computing Qij we have used reddening corrected apparent B magnitudes (where available) and 
assumed M/Lb = 20 to obtain M a jd 2 a and distances from the Tully (1988) catalog scaled to Ho — 100 
km sec _1 Mpc _1 for galaxies beyond the Local Group and the Mateo (1998) compilation for distances within 
the Local Group. For comparison with the results of RLB89 the units of the eigenvalues of Q given below 
are G x (2.32 x 10 12 M Q )/(Mpc 3 ). First, using the Local Group barycenter as origin we calculate Q for the 
effects of galaxies in the Coma-Sculptor cloud (as defined in §3.2) on the Local Group. The results are: 

Tidal Quadrupole— (Coma-Sculptor Cloud on Local Group) (n=176) 



0.0944t°;^ 3 Z = 319i|;S & = -51.3™ 



e b = -0.0122±g;°ii| I = 3W± 3 5 f 7 b = 38.2™ (16) 



-0.0822t°S Z = 43.3±i:i| b = 4.60±i~ 



Comparing these results with RLB89 we find excellent agreement both in relative amplitude, sign, and 
orientation of the eigenvalues. Wc note that the positive eigenvalue (tidal stretch) is in a direction in good 
agreement with our calculated long axis of the spatial spatial distribution of Local Group galaxies in §4.1 
while the directions of the negative eigenvalues (tidal compression) are easily identified with the previously 
determined short and intermediate axes. Thus these tides are acting to maintain the observed Local Group 
spatial distribution. 

If we now calculate the tidal effects at the Galactic center due to the Coma-Sculptor cloud and Local 
Group galaxies except the Galaxy and its companions we find 



Tidal Quadrupole— (Surrounding galaxies on the Milky Way) (n=190) 
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While the directions are similar to those in (1), (2), and (3), the amplitudes and signs are such that the 
largest positive eigenvalue (tidal stretch) is in the direction of M 31 which generally corresponds with the 
shortest spatial axis in §2.2. Similarly we find tidal compression indicated along our longest spatial axes. 
These effects would disrupt the observed satellite distribution if the tidal forces dominate. However, with 
M-B,Galaxy= —20.5 (Tully 1988) and for consistency, the above assumed M/Lb = 20, our Galaxy would 
have a mass of 0.46 x 10 12 M Q . Assuming a radius of 0.2 Mpc, the value of GM/R 3 in the above units for 
our Galaxy is then ~ 17 times the largest eigenvalue in (17), and we conclude that at least for now the 
Galaxy is dense enough to withstand the tidal forces exerted on it by its neighbors. 
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5. Discussion 

Our working hypothesis is that the spatial distribution of the companion galaxies of the Milky Way 
defines a real triaxial structure. As a measure of the degree of support for this hypothesis we identify and 
compute the difference in the direction of the long axis of the possibly related distributions from a given 
reference value. For this reference we use the long axis of the Coma-Sculptor spatial distribution in solution 
(14) and give the difference for each solution- a) for the spatial distribution of Galaxy satellites (solutions 
(1), (2), and (3)) these differences are 3.6° ±9.4, 14° ± 12, and 12° ± 12 respectively, b) the angular variation 
of radial-velocity dispersion of Galaxy satellites ((5), (6), and (7)) give differences of 2.1° ± 12, 10° ± 13, and 
21° ± 2.3. c) for the M 31 sub-group and M 31 and Galaxy satellites combined ((9) and (10)) the differences 
are 20° ± 5, and 13° ± 9.1. d) for the isolated dwarfs of the Local Group (solution (11)) the difference is 
26° ± 10. Finally for the velocity anisotropy tensor of the high-velocity clouds, solution (13), we find a 
difference of 20° ± 5.6. In addition, we note that the shortest axes of both the Local Group distribution and 
the Coma-Sculptor cloud differ by 13° ± 10 and 11° ± 2.3 from the supergalactic pole. Tidal quadrupoles, 
(16) and (17) also exhibit close alignment with both of the above principal directions with the tidal effect 
of the Coma-Sculptor galaxies acting to maintain the Local Group structure, while the tidal effects of M 31 
and surroundings are not sufficiently strong to disrupt the present structure surrounding the Galaxy. The 
outer globular clusters do not appear to be part of this same structure although they do form a nearly 
oblate, flattened system with minor axis highly inclined to the present rotation axis. 

It is the concordance of the directions of these principal axes which leads us to conclude that the above 
interpretation of what we observe in the outer regions of the Galaxy is correct, and that it is a reflection 
of structure on much larger scales. Given the above scenario we are left with having to explain why none 
of the principal axes is aligned with the Galactic pole. Modelling of galaxy formation in CDM cosmologies 
shows that the orientation of the axis of the resulting gaseous disk may bear little relation to the original 
rotation axis of the proto-galaxy because angular momentum is transfered from the inner gaseous disk to 
the dark matter dominated outer regions during the gravitational interaction with the dark matter clumps 
as they merge together. Katz & Gunn (1991) find, for example, a deviation of ~ 30° between disk axis 
and rotation axis of the original distribution of lumps. In our case this misalignment angle is ~ 12 — 30° 
assuming that the original rotation axis is the major axis. Thus we have a picture that appears plausible 
though somewhat difficult to falsify due to the small number of satellites. However, the possibility that 
the distribution of dark matter in the Galaxy may be more like that of the outer satellites is potentially 
testable. Further, there may already be some observational precedents. In addition to anticipating the 
results found here for the Milky Way, Holmberg (1969) showed on the basis of visual inspection of Palomar 
Sky Survey plates that galaxy satellites within <~ 40 h _1 kpc of many parent spirals have a tendency to be 
located along the disk minor axis. Further, Zaritsky et al. (1997) have shown that the observed distribution 
of satellite galaxies of spiral primaries out to <~ 375 h^ 1 kpc is asymmetric and is also elongated along the 
disk minor axis. 

Given the apparent ubiquitousness of this 'Holmberg' effect combined with what we have found above, 
the following galaxy formation scenario suggests itself. Early on, a preferred axis is established in the 
direction of dominant motions (such as along a filament as above) and this axis also becomes the rotation 
axis (this is certainly not unambiguously determined here, but our results for Vq are consistent with this 
idea) . This axis also serves as a conduit along which the bulk of the protogalactic lumps are and may still 
be accreted. The Katz and Gunn (KG) mechanism then causes a misalignment between this axis and the 
axis of the gaseous disk as it forms by dissipation at the center. The KG mechanism would ensure the 
observed 'stochasticity' in the degree of misalignment. Evidence for this galaxy formation scheme should 
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be revealed in the increasingly sophisticated numerical simulations of structure formation. 

The author wishes to thank Dr. Mike Hudson for many useful discussions and suggestions on this 
subject, the anonymous referee for constructive comments, and Stephen Gwyn for help with the figure. He 
is also wishes to acknowledge financial support from an NSERC of Canada operating grant. 
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Fig. 1. — The spatial distribution of 10 Galactic satellites (closed circles) and 11 M 31 satellites translated 
to the Galactic origin (open circles) plotted in Galactic coordinates. The axes are 500 kpc in length, and 
the ellipse at the center shows the distortion due to projection of a circle of radius 25 kpc in the Galactic 
plane. The quantitative description of this distribution is given in (10). 



